Ferromagnetic spinel CuCr2Se4 studied by Raman spectroscopy and lattice dynamics 

calculations 
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The lattice dynamics of the ferromagnetic spinel CuCr2Se4 (Tc ~ 430 K) was studied experi- 
mentally by measuring the Raman spectra and theoretically by calculation of zone center phonon 
frequencies within a shell model. All Raman allowed modes (^ig + Eg + 3F2g) were identified 
and assigned to specific atomic motions. The relative intensity of the Raman lines varies strongly 
with excitation photon energy between 1.58 and 2.71 eV, but no significant phonon anomalies are 
observed near Tc. 
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I. INTRODUCTION 

In the recent years there is a revived interest in 
the properties of ferromagnetic chromium chalcogenide 
spinels MCr2X4 (M=divalent element, X=S,Se,Te). 
These materils are prospective candidates for spin-based 
electronic (spintronic) applications as the strong interac- 
tion between the electronic and spin subsystems results 
in drastic changes of the electronic transport and optical 
properties near the Curie temperature, Tq.^^^ CuCr2Se4 
{Tc = 430 K) attracts particular attention as it is fer- 
romagnetic at room temperature and some theoretical 
models predict that with suitable doping it may become 
half-metallic - being an excellent metal for one spin chan- 
nel and excellent insulator for the other spin channeli^ 
While the structural, magnetic and electric properties of 
CuCr2Se4 have been subject of numerous studies, there 
are to our knowledge no reports on its lattice dynam- 
ics and possible effects of spin-phonon coupling. In this 
work we present and analyze the experimental Raman 
spectra of CuCr2Se4 single crystals in close comparison 
with predictions of lattice dynamics calculations within 
a shell modes. The Raman mode intensities exhibit res- 
onant behavior for excitation photon energies between 
1.58 and 2.71 eV, but no significant anomalies in phonon 
parameters are observed near the magnetic transition at 
Te- 



ll. SAMPLES AND METHODS 



after which the tube was cooled down and broken to col- 
lect the single platelet-like crystals with (111) surfaces 
(see Fig.l) deposited in the cooler zone. 

The Raman spectra were measured with five laser lines 
(784, 633, 515, 488, and 458 nm) using a T64000 spec- 
trometer equipped with microscope and heating stage. 
In temperature-dependent studies the laser power at the 
probe spot (3-4 /um diameter) was kept below 0.1 mW to 
assure negligible local laser heating. 




FIG. 1: (Color online) CuCr2Se4 surface topography. The 
laser focus spot is also seen. 



Single crystals of CuCr2Se4 were grown using the 
Chemical Vapor Transport (CVT) technique following 
a procedure similar to that described in Refjsl. Briefly, 
polycrystalline CuCr2Se4 was used as the precursor with 
iodine gas as transport agent placed in a sealed quartz 
tube (length: 15 cm; diameter: 1.5 cm). The source 
temperature was maintained at 870°C while the crystal 
growth at the other end of the tube occurred at 770° C. 
The reaction was carried out for a period of 80 hours. 



III. RESULTS AND DISCUSSION 

The structure of MCr2X4 is described by the Fd3m 
(No. 227) space group. The primitive cell contains 14 
atoms and the reduction of the 42-dimensional represen- 
tation r at fc = into irreducible representations gives^: 
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ig 



En + Fi, + 3F2„ + 2Au + 2Eu + 5Fi„ + 2^2« 



TABLE I: Polarization selection rules for backward scattering 
from a (111) surface. 
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The Aig, Eg, and the three F2g modes are Raman-active, 
Four of the five Fiu modes are IR-active and one is an 
acoustic mode. The Fig, 2A2u, 2F„, and 2F2u modes 
are silent. Simple calculations for backward scattering 
from a (111) surface (Table I) show that the scattered 
intensity should not depend on the rotation of the crystal 
surface around the propagation direction of the incident 
light, but only on the angle /3 between the incident el 
and scattered el polarizations. 



TABLE 11: Comparison of Raman mode frequencies (in cm ^) 
of CuCr2Se4, ZnCr2Se4, CdCr2Se4, and HgCr2Se4 at 300 K. 
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633 nm(1.96 eV), 515 nm(2.41 eV), 488 nm(2.54 eV), and 
458 nm(2.71 eV), normalized to the intensity of the Eg 
mode, are compared in Fig. 3. The most obvious change 
with increasing excitation photon energy is the increase 
of the relative intensity of the F2g modes near 170 and 
220 cm~^. Similar resonant behavior has been reported 
for CdCr2Se44 



A. Raman spectra 

Fig. 2 shows the Raman spectra of CuCr2Se4 ob- 
tained with 633 nm excitation at room temperature with 
parallel(ei || el,) and crossed(ei _L el) scattering configu- 
rations. Based on the selection rules of Table I the sym- 
metries of the Raman lines at 108 cm~^ (-^23), 144 cm~^ 
{Eg), 170 cm-1 (Fzg), 220 cm"! (Fjg), and 227 cm"! 
(Aig), are unambiguously determined. The Raman 
mode frequencies of CuCr2Se4 are compared to those of 
isostructural ZnCr2Se4^, CdCr2Se4«-® and HgCr2Se4— in 
Table II. 
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FIG. 3: (Color online) Raman spectra of CuCr2Se4 with 
784 nm(1.58 eV), 633 nm(1.96 eV), 515 nm(2.41 eV), 
488 nm(2.54 eV), and 458 nm(2.71 eV) excitation. 

The temperature-dependent Raman spectra between 
300 and 460 K as obtained with 633 nm excitation are 
shown in Fig. 4. The variations of the phonon parameters 
do not exhibit detectable anomalies near the magnetic 
transition that may be considered as manifestation of 
strong spin-phonon coupling. Here again the results are 
similar to those for CdCr2Se4 and HgCr2Se4i2ii^ 



FIG. 2: (Color online) Raman spectra of CuCr2Se4 with par- 
allel (HH) and crossed (HV) scattering configurations. 

The spectra obtained with 784 nm(1.58 eV), 



B. Lattice Dynamics Calculations 

As it is clear from the experimental results, the polar- 
ized Raman measurements allow complete discrimination 
of the Raman-active modes of different symmetry. In the 
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FIG. 4: (Color online) Temperature-dependent Raman spec- 
tra of CuCr2Se4 near the ferromagnetic transition (Tc = 
430 K). The inset shows the variations with T of the Aig 
frequency. 



TABLE III: Symmetry-adapted displacement vectors of the 
Se(u, u, u) atom in the Aig and Eg modes. The displacement 
vectors for the other atoms in the unit cell can be obtained 
by application of the corresponding symmetry transformation 
matrices.^ 



Mode symmetry 


Normalized displacement vectors 


Aig 


e = (l,l,l)/x/3 
el = (2, -1,-1)/^ 
e1 = (0,-hl,-l)/^ 
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FIG. 5: (Color online) Tetrahedral displacements of the Aig 



Eg and Aig vibrations only Se atoms are involved. The 
displacement vectors of the Se atoms in these modes are 
determined uniquely on the basis of symmetry consider- 
ations (see Table III). The corresponding displacement 
patterns are shown in Fig. 5. The Aig mode is a symmet- 
ric radial breathing of the CuSe4 tetrahedra while the Eg 
mode is a result of Se displacement tangential to the Cu- 
Se bonds. In i^2g modes, however, a strong admixture 
of Cu and Se vibrations could be expected due to the 
comparable atomic masses of copper and selenium Thus, 
lattice dynamics calculations are necessary in order to 
obtain detailed information for the atomic displacements 
in these modes. 

Lattice dynamical calculations for isostructural 
MCr2X4 chalcogenides (M=Zn,Cd,Hg; X=S,Se) re- 
ported before are based on valence force field a^°i^^ or a 
rigid- ion model (RIM). — The valence force constants are 
of limited transferability in compounds with significant 
ionic bonding. As noted in Ref. [hI, the character of the 
various phonons as relating to eigenvectors and potential 
energy distribution differs to a much larger extent than 
expected for isostructural compounds. The rigid-ion 
model proposed in Ref. [12 is based upon pairwise 
cation-anion and anion-anion short range potentials. 
The potential parameters, however, are mapped onto a 
set of effective stretching and bending force constants 
constants. The reported effective force constants do 
not display a clear trend over the series of different M 
atoms, and their transferability in the case of M=Cu 
is not evident. Moreover, RIM does not account for 
the atomic polarizabilities and is unable to reproduce 
other important lattice properties as the high-frequency 
dielectric constant, the LO-TO splitting, as well as the 
off-diagonal elastic constants. 

For the reasons outlined above we performed lattice 
dynamics calculations for CuCr2Se4 in the frame of the 



shell model (SM), which is the minimal extension over 
RIM capable to account for the ionic polarizability and 
to predict correct LO-TO splitting. In order to reduce 
the number of adjustable model parameters some phys- 
ical simplifications were adopted. First, a valence shell 
was considered for the Se atoms only while Cu and Cr 
were treated as rigid ions. Second, the van der Waals at- 
tractive interaction was considered to act only between 
Se shells, while it was neglected for the Cu (Cr) core - Se 
shell pairs. These assumptions are justified by the much 
higher polarizability of the Se~^ ion compared to that 
of the transition-metal ions. The rigid-ion approxima- 
tion for transition-metals is a common approximation in 
the shell-model calculations on transition-metal oxides J^ 
The short-range interatomic interactions were modelled 
by a Buckingham potential: V{r) — Aexp{—r/p) — C/r^, 
where a non-zero van der Waals constant C was retained 
for the Se shell - Se shell pairs only. 

Even with the above simplifications the number of ad- 
justable parameters exceeds the number of observables 
for CuCr2Se4, the lattice constant a, the Se fractional 
coordinate u and the five Raman-mode frequencies. For 
this reason, the A, p and C parameters for the Se-Se in- 
teraction, the Se core and shell charges X and Y, and the 
Se core-shell spring constant k were adjusted initially by 
fitting of the lattice constant, static and optical dielectric 
constants, the LO and TO phonon frequencies, and the 
elastic modulus in ZnSe. As a next step these parameters 
were transferred to CuCr2Se4 and the parameters of the 
Cu-Se and Cr-Se core-shell interactions were optimized 
in order to fit the lattice parameter, the Se fractional 
coordinate u and the frequencies of the Raman-active 
modes. The list of as obtained SM parameters is shown 
in Table IV. 

A comparison between the calculated and the experi- 



TABLE IV: Adjusted shell- model parameters for CuCr2Se4. 



Atom 


Core charge 
X 


Shell charge 
Y 


Core-shell 
spring constant 

k (eV/A^) 


Cu 
Cr 
Se 


1.00 
2.30 
2.19 


-3.59 


69.7 


Atomic pair 


A{oV) 


P(A) 


C (eVxA«) 


Cu core - Se shell 
Cr core - Se shell 
Se shell - Se shell 


6827.4 
13817.5 
1482.2 


0.2364 
0.2513 
0.3440 






136.3 



TABLE V: Cell parameters, frequencies (in cm ^) and atomic 
motions. 
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FIG. 6: (Color online) Atomic displacement of the Raman 
modes of CuCr2Se4. 



mental lattice parameters and Raman-mode frequencies 
is shown in Table V. The calculated atomic displace- 
ments for the three Fig modes are depicted in Fig. 6. The 
lowest-frequency mode at 108 cm"-"^ consists of an almost 
pure Cu translation within a practically static Se cage. 
This assignment corroborates the fact that the mode fre- 
quency scales with the mass M of the A atom (A — 
Cu, Cd, Hg) as Af^i/2 (see Table-II). As an additional 
piece of evidence, we performed phonon calculations by 
substituting Cu with Cd and Hg and keeping the same 
values for the interatomic potentials. We obtained the 
frequencies of the lowest F^g mode at 83 cm~^ for A = 
Cd and 63 cm~^ for A — Hg in perfect agreement with 
the experimental values. The medium-frequency mode 
at 173 cm~^ can be described as Cu translation against 
the Se cage, which moves in the opposite direction. The 
highest-frequency vibration at 200 cm"^ is a superposi- 
tion of an in-phase stretching of the Cu-Se bonds and 



C. Resonant Behavior of Raman Scattering 

The Raman scattering cross-section S{huj) in non- 
transparent materials is related to the combined elec- 
tronic density of states. As a rule S{hLo) is proportional 

^^ ^^^ I dhLj \'^ '^^'^ ^^^ maxima near interband gaps. 
The experimentally observed dependence of the Raman 
intensity on the excitation photon energy I{fvjj), however, 
is additionally modified by the spectral properties of the 
Raman setup and the absorption and reflection losses, 
which occur for the incident and scattered light. The 
detail analysis of resonant Raman scattering in closely 
related CdCr2Se4 has shown that S{huj) has a maximum 
near 2 eV<^ It has also been shown in that study that the 
intensity of the Fig Raman lines, if normalized to the Eg 
intensity, strongly increases with photon energy between 
1.5 and 2.8 eV thus demonstrating the dependence of 
scattering cross section on the phonon symmetry. As it 
follows from Fig. 3, similar enhancement of the relative in- 
tensity of the Fig modes is observed for CuCr2Se4. Fig. 7 
shows in more detail the variation with huj) of the nor- 
malized intensity for the 170 cin~^ (Fig), 220 cm~^(i^2g), 
and 227 cin~^{Aig) modes. 

For tentative explanation of the resonant enhancement 
of the Fig intensity let us consider the ab initio electronic 
structure calculationsi^ii^ for CuCr2Se4. Following these 
calculation the electronic density of states (DOS) above 
the Fermi level is dominated by the Cr d states, while 
those below the Fermi level by the Cu d and the Se 
p states. An optical transition between Cu and Cr d 
bands can be disregarded due to the large spatial separa- 
tion between the Cu and Cr ions. Therefore, a possible 
candidate for an optical transition could be an electron 
transfer between the Se p and Cr d orbitals. Transition 
energies of 3.2 eV between Se pi/i and Cr d^/i states, 
and 2.84 eV between Se p^/i and Cr d^/i states could 
be expected on the basis of calculated centers of gravity 
of the corresponding bands^^. The latter value is close 



the experimentally observed maximum at 2.75 eV of the 
optical conductivity of CuCr2Se4 reported by Ongushi et 
aLifi 



sentations for the 16d (Cr) and 32e (Se) positions gives: 
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r(Cr) = V4lg + F2g 



(1) 
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Therefore, among Raman-active modes only Ai^ and F2g 
modulate the Se and Cr orbital energies and respectively 
the energy of the p ^ d transitions. This fact can be un- 
derstood by inspecting the atomic displacement patterns 
in Fig. 6. In the Aig and F2g modes the Se displacements 
have identical radial components relative to the Cr atom 
and, therefore, give identical variations of the Madelung 
potential at the Cr site. For the Eg mode, however, the 
radial components of the Se displacements are of alter- 
nating signs leading to cancelation of the Madelung po- 
tential variations at the Cr site. By following our lattice 
dynamics calculations it is easy to note that the reso- 
nant enhancement of the ¥2g modes increases with the 
increase of the radial character of the vibration. Inter- 
estingly, the Aig mode shows an enhanced normalized 
intensity also with 784 nm(1.58 eV) excitation, which is 
in the energy range of the Cr d — d transitions.-—. 



FIG. 7: (Color online) Variation of the normalized Ra- 
man intensity of the 170 cm~ {F2g), 220 cm~ {F2g), and 
227 cm~^(j4ig) modes with excitation photon energy. 

The lattice dynamics calculation can help us under- 
stand why Aig and i^2g niodes exhibit stronger coupling 
to Se p - Cr d transitions than the Eg. The coupling of a 
Raman-active phonon to an optical transition is mainly 
due to two mechanisms: (i) modulation of the dipolar 
transition moment, and (ii) modulation of the transition 
energy by the atomic displacements. As a rule the mod- 
ulation of the transition energy has stronger resonant be- 
havior and gives a dominant contribution to the Raman 
scattering intensity under resonant conditions. We will 
consider the optical p — ?> d transitions in a simplified 
tight-binding picture assuming that the transition energy 
is roughly equal to the difference between the on-site en- 
ergies of the Se p and Cr d orbitals. Since the on-site 
energy is a scalar quantity, it transforms according to 
the permutational representation of the given atomic po- 
sition. The decomposition of the permutational repre- 



IV. CONCLUSIONS 

The Raman spectra of the ferromagnetic spinel 
CuCr2Se4 {Tc = 430 K) was studied experimentally and 
analyzed theoretically by calculation of lattice dynam- 
ics. The normalized Raman scattering intensity for Aig 
and F2g exhibits resonant behavior similar to that known 
for the isostructural CdCr2Se4, but no significant phonon 
anomalies are observed near Tc- 
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